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We introduce the concept of controlling the nonlinear response of the metamaterial by altering its internal 
structure. We experimentally demonstrate that we can tune the nonlinear frequency shift for two coupled 
split-ring resonators by changing their mutual position. This effect is achieved through modification of the 
structure of the coupled resonant modes and their interaction with the incident field. We explain these results 
with the aid of numerical modeling, which shows that by offsetting the resonators we control the maximum 
currents through the nonlinear driving elements which affects the nonlinear response of the system. 
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The concept of engineering the electromagnetic prop- 
erties of a material to possess a negative refractive index 
was first put forward by Veselago in 196P. The first 
experimental confirmation of this idea was accomplished 
by Smith et al. in 200G 2 , using a metamaterial consist- 
ing of a combination of split-ring resonators (SRRs) to 
provide a negative permeability^, and thin metal wires to 
provide a negative permittivity. More recently, metama- 
terials composed of only SRRs have been used to produce 
optical activity^, offering an alternate route to negative 
refractiorP. 

The effective properties of these metamaterials depend 
both on the properties of the individual SRRs, and on 
the arrangement of the SRRs in the lattice, as strong 
interactions occur between the rings in a lattice due to 
their complex near-field patterns. This allows us to sig- 
nificantly control the properties of the metamaterial by 
changing the internal structure of the composite mate- 
riaP. In order to understand the resulting responses of 
the metamaterial, we need to understand the interactions 
between the individual SRRs. This is a major motivation 
behind the studies involving the relative positions of two 
SRRs, including rotating therrpl, and shifting them with 
respect to each other^. 

The linear properties of individual metamaterial ele- 
ments can also be controlled by an external signal, by in- 
serting nonlinear inclusions such as semiconductors^ or 
diodes^ into the rings. Many tunable inclusions can also 
respond to the high-frequency incident wave, thus they 
form the basis of nonlinear metamaterials, whic h typ- 
ically have power-dependent resonant frequencies^"^. 
The nonlinear response of metamaterials can be much 
stronger than that of natural materials, due to their res- 
onant nature, and the local field enhancement which oc- 
curs at certain "hot spots" . 

The strength of the local field depends on the design 
of the resonator, and the choice of an optimal location 
to place the nonlinear element. However it also depends 
strongly on the coupling of the wave to the external field, 
which along with losses determines the quality factor of 
the resonator. Since modifications of the lattice param- 
eters also influence the quality factor of the resonances® , 
they should affect the local-field enhancement. In this 
Letter we experimentally demonstrate, for the first time 
to our knowledge, a control of the nonlinear response of 
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FIG. 1. (a) Schematic of the lateral shifting of SRRs with 
inserted diodes; Sa is a shift between two rings, (b) The 
measured absorption curves for 5a — 3.75 mm (blue solid) and 
7.5mm (green dashed), at -20dBm input power. Symmetric 
(ujs) and antisymmetric (ojas) modes are highlighted, (c) 
Resonant frequencies of symmetric (black circle/solid) and 
antisymmetric (red cross/dashed) resonances, as determined 
both experimentally (markers) and numerically (lines) for the 
linear case. 



two broadside-coupled split-ring resonators by modifying 
the offset between them, and explain this response by 
studying how shifting the resonances changes the maxi- 
mum currents in the resonators. 

First, we perform microwave experiments with a pair of 
SRRs, having an offset Sa between their centers, as shown 
schematically in Fig. [fja). The copper rings are 3.75 mm 
in outer radius, 3.25 mm in inner radius, and have a gap 
width of I mm. These rings are printed on the opposite 
sides of 1.6 mm thick FR4 circuit board. Each ring has a 
second gap of 0.4 mm opposite the first gap, across which 
a Skyworks SMV1405-079 series diode is soldered. The 
samples were placed inside WR-229 rectangular waveg- 
uide, oriented such that the incoming magnetic field is 
perpendicular to the loops. 

The transmission and reflection from the sample are 
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measured using a Rohde and Schwarz ZVB-20 vector net- 
work analyzer. We measure the absorption spectrum, 
which describes the frequency of the maximum current 
distribution in the rings, given by 1 — |52i| 2 — |Sii| 2 , 
where S21 and Sn are the transmission and reflection 
coefficients, respectively. The absorption curves are mea- 
sured experimentally for the values of 5a between and 
7.5 mm, in 0.375 mm increments, at powers -20dBm and 
f 5 dBm. Numerical results are calculated using CST Mi- 
crowave Studio. Each diode is replaced by a series RLC 
circuit with 2.67 pF capacitor, 0.851 resistor and 0.7 nH 
inductor. Such a circuit describes the diode in the linear 
regime^. 

Figure [TJb) shows the measured absorption curves for 
5a equal to 3.75 mm, and 7.5 mm, at -20 dBm. There 
are seen clearly two resonant modes, symmetric and an- 
tisymmetric, which are defined by the relative directions 
of the currents on the rings. The se modes are both 
tunable through a change of 8a®®, which is evident, 
as when 5a = 3.75 mm the symmetric (antisymmetric) 
mode is the lower (higher) resonant mode, while when 
da = 7.5 mm it has increased (decreased) to become the 
higher (lower) resonant mode. The frequencies of max- 
imum absorption for both modes are found and plotted 
as a function of 5a in Fig.[T|c) (markers), along with the 
numerical equivalents (solid lines). We observe excellent 
agreement between the experimental and numerical re- 
sults. Consistent with the findings in Ref. HI there are 
strong responses in both the modes, with the resonant 
frequency of the symmetric mode increasing with an in- 
crease in 5a. 

As seen in Fig. [l|b) , the symmetric mode is much more 
strongly excited than the antisymmetric mode. In fact, 
the antisymmetric mode is not visible for low values of 
5a, as seen in Fig. |T](c), as the coupling to the incident 
wave is very weak. However, we have enough points to 
see that the antisymmetric mode is still tunable via the 
lateral shifting. 

In a single SRR, the resonant frequency is determined 
by u> = 1/v 7 LC , where L is inductance, and C is ca- 
pacitance. Introducing the diode adds a nonlinear ca- 
pacitance Cy in series with the capacitance of the SRR. 
By adding the second SRR, we introduce coupling be- 
tween the rings, which causes a split in the resonant fre- 
quency, resulting in the two modes previously discussed, 
ojs and (±>as- By changing the intensity of the incom- 
ing microwaves, we are causing a nonlinear change in 
the capacitance of the diodes, which changes the effec- 
tive capacitance of the system, and shifts the resonant 
frequencies. This results in having an added degree of 
freedom for manipulating our system. 

In our structure, the nonlinear response is quantified 
by a shift of the resonant frequency with increasing input 
power from -20dBm to 15dBm. The relative difference 
between the low- and high-power resonant frequencies 
for each value of 5a is shown in Fig. [2j The inset in 
Fig. [2] shows the absorption of the symmetric peak for 
both powers when 5a = 3.75 mm. We observe a signif- 
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FIG. 2. Relative nonlinear shift in symmetric (black circles) 
and antisymmetric (red crosses) resonant frequencies for shifts 
(5a) . The dashed lines are second order polynomial fits to the 
data. The inset shows the nonlinear shift of the absorption 
at the symmetric mode when 5a = 3.75 mm. The solid blue 
line is at -20dBm, and the green dashed line at 15dBm. 



icant change of the nonlinear frequency shift as a func- 
tion of the linear offset. It is clear that the nonlinearity 
is much stronger for the symmetric resonance when the 
rings are closer together. It can be seen that the anti- 
symmetric mode is also tuned by the increasing power, 
though due to the excitation we use, we cannot excite 
this mode for smaller offsets. 

We can explain the nonlinear tuning by looking both 
experimentally and numerically at some of the linear as- 
pects of the system. As the nonlinear shift is due to the 
currents excited in the rings (which also leads to absorp- 
tion) , it would be expected the higher the absorption, the 
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FIG. 3. (a) Experimental (markers) and numerical (lines) 
maximum absorptions for both symmetric (black circle/solid) 
and antisymmetric (red cross/dashed) modes, (b) Numer- 
ically calculated absorption coefficients at the resonant fre- 
quencies for the case of a lossless substrate. 
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FIG. 4. Experimental (markers) and numerical (solid line) Q 
factors of the symmetric mode. 



FIG. 5. (a) Calculated voltage over the capacitive compo- 
nent of the diode, (b) Resulting capacitance (presented as 
l/sqrt(C)). 



larger the nonlinear response. Therefore the maximum 
absorption for each value of 5a is plotted in Fig. [3^a) 
for both modes, showing good agreement between exper- 
imental and numerical results. However, by comparing 
this figure with Fig. [2j we see increasing absorption with 
decreasing nonlinear shift. This effect appears due to 
the fact that, offsetting the rings, we not only modify 
the current amplitudes in the rings, but also change the 
field distribution in the substrate. As a result, the total 
absorption presented in Fig. 3]ja) contains contribution 
from losses in both metal and in dielectric. At the same 
time, the nonlinear response in our case is caused by 
currents in metal only. The losses within the FR4 dielec- 
tric circuit board increase with 5a due to increased field 
confinement within the board, and change in these losses 
dominate changes of the absorption in metal. In Fig.^b) 
we show the results of a simulation which neglects dielec- 
tric losses. We can see that these maxima now decrease 
as 5a increases, agreeing with the trend in the nonlinear 
response. The difference in trends between Figs.[3]ja & b) 
confirms explained above effect of dielectric losses. 

The amplitude of the currents in the system depends 
on the Q factor of the resonator. Because of this we would 
expect to see some correlation between the changing non- 
linear shift, and the Q factor. By comparing Fig. [2] with 
Fig. |3J where we have plotted the Q factor as a function 
of 5a for the symmetric mode, we can see that this is in- 
deed the case, and decreasing quality factor leads to the 
decrease of the nonlinear response. 

The physics of the nonlinear shift of the resonant fre- 
quency was studied in Rcf. 12 -", where it was shown that 
the rectified AC voltage across the varactor provides self- 
biasing mechanism. By increasing the amplitude of the 
electromagnetic wave, we increase alternating voltages 
in the SRR, generating a larger DC bias voltage, which 
reduces the varactor capacitance. To study this effect, 
we performed numerical simulations, and monitored the 
voltages on each RLC circuit representing a varactor. As 
a larger voltage leads to a stronger decrease in the varac- 
tor capacitance, we would expect to see a similar trend 
between the voltages recorded at the resonant frequency, 
and the nonlinear shift. We want to note that the voltage 
across the capacitance, and not the whole RLC circuit 



should be calculated. Corresponding DC voltage at the 
resonant frequency is plotted for both resonant modes in 
Fig. Ha). 

By comparing Figs. |2j ^b), and[5ja), we can see an 
overall correlation, where as 5a increases there is a de- 
creasing trend in the nonlinear response, the maximum 
lossless absorptions, and the calculated voltages, for the 
symmetric mode. This enables us to explain the nonlin- 
ear response of the system in a meaningful way. 

As the nonlinear response should be directly propor- 
tional to the inverse of the square root of the capacitance, 
we then calculated this capacitance using the known volt- 
age across the diode. From the data sheets for the diode, 
we can obtain the capacitance (Cy) for a given applied 
voltage (Vr) using the equation^ 
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where Cjo is the zero bias junction capacitance, Vj is the 
junction potential, M is the grading coefficient, and Cp 
is the package capacitance. For the specific diode series 
we have used, these values are Cjq = 2.92pF, Vj = 0.68V, 
M = 0.41, and C P = 0.05pF. The inverse of the square 
root of the resulting capacitance is shown in Fig. [5jb) , as 
a function of 5a. As expected, this plot closely resembles 
the plots in Fig. [5ja) and Fig. [2j The dependence of 
the diode capacitance on 5a changes the response of the 
system to linear tuning of the resonant frequencies, since, 
as discussed earlier, the resonant frequency is dependent 
on the combined capacitance of the diode, and the SRR. 

When the two resonances overlap, it is difficult to ac- 
curately distinguish them in the data. Therefore, the 
experimental values for 5a in the shaded regions have 
been excluded for Figs.[ljc) -[3^a). 

In conclusion, we have demonstrated that by shifting 
two coupled split-ring resonators relative to each other, 
we can significantly control the nonlinear properties of 
both the symmetric and antisymmetric resonant modes 
in the system. We have also found that the tuning of 
the nonlinear response can be explained and predicted 
numerically by studying the voltages generated across the 
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diodes. We expect that our results will stimulate further 
work in controlling and designing nonlinear properties of 
metamaterials. 
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